Law for the distance between successive earthquakes 
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Abstract. Discovery of a new law for the three-dimensional spatial distance between 
the foci of successive earthquakes is reported. Analyzing the seismic data taken between 
1984 and 2001 in southern CaUfomian, it is found that the cumulative distribution of the 
distances follows the modified Zipf-Mandelbrot law, showing complexity of geometry 
of the events. 
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1. Introduction 



Though seismicity is characterized by extraordinarily complex phenomenology, some 
of the known empirical laws are remarkably simple. Celebrated examples are the Omori 
law [Omori, 1894] for temporal distribution of aftershocks and the Gutenberg-Richter 
law [Gutenberg and Richter, 1944] for the relationship between frequency and 
magnitude. These classical examples exhibit the scale-free nature of earthquake 
phenomena, and this is why seismology has been attracting much attention of 
researchers of complex systems [Bak and Tang, 1989; Olami et ah, 1992; Huang et ah, 
2000; Baket a/., 2002]. 

In a recent work [Abe and Suzuki, 2002], we have studied temporal complexity of 
earthquakes in southern California. We have analyzed the waiting time distribution 
between successive earthquakes, and have found that the seismic time series undergoes 
over a series of transitions between quasi-equilibrium states, each of which obeys the 
characteristic distribution termed the "^-exponential distribution". The ^-exponential 
distribution is a generalization of the Zipf -Mandelbrot distribution [Mandelbrot, 1983], 
and is frequently encountered in studies of complex systems at their quasi-equilibrium 
states. It can be derived, via maximum entropy principle, from Tsallis's generalized 
entropy [Tsallis, 1988; Tsallis et ah, 1998; Abe, 2002]. The associated generalized 
statistical theory, referred to as nonextensive statistical mechanics [Abe and Okamoto, 
2001; Kaniadakis et al., 2002; Gell-Mann and Tsallis, in press], is known to be in 
conformity with the principles of macroscopic thermodynamics, and is currently under 
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active investigation. 

Regarding spatial complexity in seismicity, it may be worth pointing out [Okubo and 
Aki, 1987; Marsan et al., 2000] that the distributions of the faults and the foci have the 
fractal structures. It is of fundamental interest and importance to predict where the focus 
of the next shock will occur after earlier shocks. Science does not seem to be able to 
make such a prediction yet but, as a step toward it, it may be worth clarifying the 
statistical property of the geometric distances between successive earthquakes. 

In this paper, we report the discovery of a new empirical law for the nature of spatial 
geometry of earthquakes. In particular, we analyze the statistical property of the three- 
dimensional distances between successive earthquakes. We show that the resulting 
cumulative distribution of the distances is described by the ^-exponential distributions. 
This result exhibits spatial complexity of earthquakes in a novel manner. 

2. Modified Zipf -Mandelbrot Law and ^-Exponential Distribution 

In this section, we discuss the statistical foundations for the Zipf -Mandelbrot law and 
its modification in view of nonextensive statistical mechanics. 

Nonextensive statistical mechanics is a scale-invariant theory, which can explain the 
statistical properties of a variety of complex systems at their quasi-equilibrium states. It 
can be formulated by the maximum entropy principle. Previously, in seismology, 
maximum entropy principle was examined for describing, for example, the frequency- 
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magnitude distribution [Main and Burton, 1984], in which the ordinary Boltzmann- 
Gibbs -Shannon entropy was employed. Here, the entropy functional to be considered is 
a generalized one termed the Tsallis entropy. 

Suppose the three-dimensional distance, r, between the foci of successive 
earthquakes be a fundamental random variable to be measured. p(r)dr stands for the 
probability of finding the value of the distance in the range, [r, r + dr]. Therefore, the 
normalization condition reads 

]drp(r) = L (1) 



In reality, p (r) has a compact support with the finite maximum value, r^^^ . 
The Tsallis entropy is given as a functional of p(r) as follows: 



1-q 



Vo 



] — [op(r)r-l 



(2) 



where q and (T are the positive entropic index and a scale factor of the dimension of 
length, respectively. In the limit q^l, this quantity converges to the Boltzmann- 
Gibbs-Shannon entropy: 

limS^ip] = -Idr p(r)ln[(7p(r)]. (3) 
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Henceforth, cr is set equal to unity for simplicity. The Tsallis entropy is known to share 
a lot of common properties with the Boltzmann-Gibbs-Shannon entropy, including 
concavity and fulfillment of the //-theorem. However, additivity is violated. For a 
system composed of two statistically independent systems, A and B, the Tsallis entropy 
satisfies 



which is referred to as pseudoadditivity. The last term on the right-hand side brings the 
origin of nonextensivity of the resulting generalized statistical mechanics. Therefore, the 
value, indicates the degree of nonextensivity in a complex system under 

consideration. 

In nonextensive statistical mechanics [Tsallis et al., 1998; Abe and Okamoto, 2001], 
the definition of the expectation value should also be generalized to the "^-expectation 
value", according to the generalization of the entropy. The expectation value of the 
distance is given by 



S^(A, B) = S^(A) + S^(B) + q)S^(A)S^(B), 



(4) 




(5) 







where P^ir) is the escort distribution [Beck and Schldgl, 1993] defined by 
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\dr'[p(r')Y 

J 



(6) 



Now, maximization of the Tsallis entropy under the constraints on the normalization 
in Eq. (1) and the ^-expectation value of the distance in Eq. (5) leads to the following 
optimal distribution [Tsallis et al, 1998; Abe and Okamoto, 2001]: 



--(r-<r> ) 



(7) 



Z^=\dre^ 





-^{r-<r> ) 
\ c " J 



(8) 



c = ]dr{p,{r)V 



(9) 



Here, P is the Lagrange multiplier associated with the constraint on the ^-expectation 
value of the distance. e^{x) denotes the "^-exponential function" defined by 



[\ + {\-q)x\ "^'-'''> (l + (l-q)x>0) 
(l + (l-q)x<0) 



(10) 



whose inverse is the "^-logarithmic function" 



6 



ln,(x) = -^(^^-^-l). (11) 
l-a 



In the limit q^l, the ^-exponential and ^-logarithmic functions tend to the ordinary 
exponential and logarithmic functions, respectively. 

Pair) in Eq. (7) is a power-law distribution unless in the limit ^ — > 1. It has a long 
tail if q>l, whereas there appears the cut-off at = / (I- q) if 0<q<l, where 
rQ=[c + {l-q)P<r>^]/ P, which is positive even if q>l [Abe and Rajagopal, 
2000a]. Therefore, r^^^ = if 0<q<l, whereas r^^ can be arbitrarily large if q>l. 
Using this r^, Poir) is rewritten as follows: 



Po(r) = -z^ (12) 

I dr'e^(-r'/r,) 



This is referred to as the "^-exponential distribution". It is a generalization of the Zipf- 
Mandelbrot distribution. (The standard Zipf -Mandelbrot distribution corresponds to the 
case when ^ > 1.) The exponential distribution is seen to be recovered in the limit 
q^l. 

A point of crucial importance in nonextensive statistical mechanics is that the 
quantity to be compared with the observed distribution is not Poir) but its associated 
escort distribution [Abe and Rajagopal, 2000b]. Accordingly, the cumulative 
distribution should be defined by 
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P(> r) = j Jr' P^ir'). 



(13) 



Using Eqs. (6) and (12), we obtain 



P(> r) = e^(-r/ r,). 



(14) 



In what follows, we shall show by data analysis that the statistical property of the 
three-dimensional distances between the foci of successive earthquakes is described 
extremely well by the distribution in Eq. (14) with 0<q<l, exhibiting spatial 
complexity of earthquakes. 

3. Distribution in Southern California 

In this section, we present the data analysis and its comparison with the modified 
Zipf -Mandelbrot law. 

We have analyzed the earthquake catalog made available by the Southern CaUfomia 
Earthquake Data Center (http://www.scecdc.scec.org/catalogs.html) covering the period 
between 00:25:8.58 on 1 January 1984 and 23:44:2.81 on 31 December 2001 in the 
region spanning 29°15.25'N- 38°49.02'N latitude and 113°09.00' W-122°23.55' W 
longitude. (We have taken this period since the data in 1983 are partially missing for a 
few months.) The number of events is 364867. The data contain, not only significant 
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earthquakes but also extremely weak earthquakes like magnitude 0.0. We have 
calculated the three-dimensional distances between successive earthquakes. 

In Figures la and lb, we present the log-log and semi- log plots of the cumulative 
distribution associated with the statistical frequency of the distance, respectively. The 
dots represent the observed distribution, whereas the solid line corresponds to the Zipf- 
Mandelbrot law described by the ^-exponential function. We also present the semi-^-log 
plot in Figure Ic. It is clearly appreciated that quality of the fitting is extremely high. 

4. Distribution in Japan 

For comparison with the result obtained in the previous section, we here present the 
analysis of the catalog of earthquakes in Japan made available by the Japan University 
Network Earthquake Catalog (http://kea.eri.u-tokyo.ac.jp/CATALOG/junec/monthly. 
html) covering the period between 01:14:57.63 on 1 January 1993 and 20:54:38.95 on 
31 December 1998 in the region spanning 25.851° N- 47.831° N latitude and 
126.433° E- 148.000° E longitude. The number of events is 123390. (We have limited 
ourselves to this period since before 1993 the number of the observed data per year 
turned out to be about half of the latter period. An essential difference of the catalog of 
the Japan University Network Earthquake Catalog from that of the Southern California 
Earthquake Data Center is that the former does not contain earthquakes with magnitudes 
smaller than 2, unfortunately.) 
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In Figures 2a and 2b, we present the log-log and semi-log plots of the cumulative 
distribution associated with the statistical frequency of the distance, respectively. The 
dots represent the observed distribution, whereas the solid line corresponds to the Zipf- 
Mandelbrot law described by the ^-exponential function. We also present the semi-^-log 
plot in Figure 2c. As in the case in southern CaUfomia, the fitting is seen to be very 
well. 

5. Conclusions 

We have studied the statistical property of the three-dimensional distances between 
successive earthquakes and have discovered that it obeys the modified Zipf-Mandelbrot 
law characterized by the ^-exponential distributions with q = q, less than unity 
= 0.75 ~ 0.77). This result exhibits complex spatial geometry of earthquake 
phenomenon in a novel manner. 

In our previous work on the time intervals between successive earthquakes [Abe and 
Suzuki, 2002], the associated waiting time distribution was shown to be given also by 
the ^-exponential distribution with q = q,>l (typically, q, = 1.2 ~ 1.3). It is of interest 
to notice that the duality relation, q^ + q,~l, might hold, though it is still 
hypothetical. 
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Figure Captions 



Figure la Log-log plot of the cumulative distribution of the three-dimensional 
distances between successive earthquakes in southern CaUfomia. The 
values of the index, q, and the parameter, r^, are q = 0.113 and 
= 1.79 X 10^ km, respectively. 

Figure lb Semi-log plot of the data in Figure la. 

Figure Ic Plot of the data in Figure la on semi-^-log scale. The straight line describes 
the ^-exponential distribution. The value of correlation coefficient is 
p= -0.9993. 

Figure 2a Log-log plot of the cumulative distribution of the three-dimensional 

distances between successive earthquakes in Japan. The values of the index, 
^, and the parameter, Tq, are ^ = 0.747, = 5.95x10^ km, respectively. 

Figure 2b Semi-log plot of the data in Figure 2a. 

Figure 2c Plot of the data in Figure 2a on semi-^-log scale. The straight line describes 
the ^-exponential distribution. The value of correlation coefficient is 
p= -0.9990. 
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Figure la 




Figure lb 
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